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Tij- ■ In the present paper few steps are undertaken towards the descrip- 

\^ . tion of the "qubit-qutrit" pair - quantum bipartite system composed 

I of two and three level subsystems. The computational difficulties with 

the construction of the "/oca/ unitary polynomial invariants" are dis- 
cussed. Calculations of the Molien functions and Poincare series for 
the qubit-qubit and qubit-qutrit local unitary invariants are outlined 
^ ' and compared with the known results. The requirement of positive 

■ semi-definiteness of the density operator is formulated explicitly as a 

set of inequalities in five Casimir invariants of the algebra su(6). 
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1 Introduction 



The present article discusses several computational aspects of a pure quantum 
effects in composite systems playing an important role in the modern theory 
of quantum computing and quantum information [HE]- 

The cornerstone of these latest trends is an extraordinary quantum phe- 
nomenon - the "entanglement" of quantum states. Basically, under the en- 
tanglement it is assumed an exposition of diverse non-local correlations in a 
composite multipartite quantum system, which have no classical analogue. 
From the mathematical standpoint of view characteristics of entanglement 
can be understood within the classical theory of invariants (cf. [SI II])- The 
central object in these studies is the ring of G-invariant polynomials, called 
local invariants, in elements of the density matrix with the group G con- 
sisting from the so-called local unitary transformations acting separately on 
every part of the multipartite composite system. The program of description 
of this ring for multipartite mixed states was outlined in [S] and during the 
last decade has been intensively developed. Over this time many interesting 
physical and pure mathematical results have been obtained. Particularly, for 
the simplest bipartite system of two qubits, the structure of the correspond- 
ing ring has been clarified (see e.g. O [71 IH])- However, comparative less is 
known for multipartite states, as well as for bipartite mixed states, composed 
from arbitrary d-level subsystems, the so-called qudits [HI [ID]- The reason 
is first of all in a big computational difficulties we are faced. Indeed, even 
dealing with 3-level subsystem, qutrit, the large number of independent ele- 
ments of the density matrix leads to the wide variety of the local polynomial 
invariants and makes non-effective the direct usage of the modern computer 
algebra software. 

Below, attempting to construct the polynomial ring of local invariants for 
qubit-qutrit pair, i.e. invariants against the action of SU(2)(8)SU(3) group, we 
got added evidence of the complexity of the problem. The known results [12] 
and our calculation of the Molien function and Poincare series show that 
the number of local invariants grows up significantly compared with the two 
qubits case. Nevertheless the derived information is very useful for the anal- 
ysis of the polynomial ring of local invariants. As a preliminary result we 
present here a set of linearly independent local invariants up to the fourth or- 
der constructed via trace operation from the non-commutative monomials in 
three elements of a special decomposition of qubit-qutrit density matrix. Us- 
ing the subset of the local invariants, consisting from the Casimir invariants 
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of the enveloping algebra it(su(6)) , the positive semi-definiteness of density 
matrix of qubit-qutrit pair is derived in the form of a system of algebraic 
inequalities. 

2 The SU(n) Casimir invariants 

Here the basic statements on the unitary symmetry in quantum mechanics 
and its role in the description of composite multipartite states is given. 

• Density operator and SU(n)-invariants • According to the 
conventional quantum theory, a complete information about a generic n- 
dimensional system is accumulated in the self-adjoint positive semi-definite 
density operator g with the unit trace, g G For a closed quantum sys- 
tem, this description is highly redundant, the equivalence relation between 
elements of , due to the invariance of observables under the adjoint action 
of SU(n) group 

iAdg)g = ggg'\ ^? G SU(n) , (2.1) 

guarantees that the physically relevant knowledge about quantum states can 
be extracted from the orbit space ^+ |SU(n) Q|. Relaxing for a moment 
condition of semi-definiteness, the density operator g can be expressed via 
the Lie algebra su(n) of SU(n) group [TT] : 

g=-In + Kig, G 5u(n) , i^ = -l. (2.2) 
n 

with some normalization factor R . Therefore the density operator can be 
decomposed over n^ — 1 basis elements, , of the Lie algebra su(n) 

n2-l 

0=J]6e., (2.3) 

i=l 

and any other operator A[g] , constructed from the density operator p, admits 
a representation in the graded power series: 

^(e) = + ) e. + i^Sf e,e, + ^A^ e,e,e, + . . . . (2.4) 

^The orbit space ^+ |SU(n) of SU(n) is defined as the set of all SU(n)-orbits, endowed 
with the quotient topology and differentiable structure, the subset of all the SU(n)-orbits 
with the same orbit-type forms a stratum of Cp+ |SU(n) . 
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According to the Poincare-Birkhoff-Witt theorem [13] the ordered monomials 



form a hnear basis of the universal enveloping algebra il(su(n)) of su(n) . 
Direct corollary of this theorem is that the symmetrized monomials of degree 
d in (12 ■4p span a linear spaces il'^(5u(n)) and the universal enveloping algebra 

oo 

il(su(n)) = il'^(su(n)) . 

(1=0 

as a linear space is isomorphic to a polynomial algebra in commutative real 
variables , i = 1, . . . , n^ — 1. 

Furthermore, according to the well-known Gelfand's theorem [13], the 
description of center, Z(su(n)) , of the enveloping algebra it(su(n)) reduces 
to the study of invariants in commutative symmetrized algebra 5'(5u(n)) , 
which is isomorphic to the algebra of invariant polynomials over su(n). The 
elements of center Z{su{n)) are in one to one correspondence with the SU(n)- 
invariant polynomials in n^ — 1 real variables, coordinates in 5u(n). More 
precisely, the element of il(su(n)) 



aG6r 

where &r is the group of permutation of 1, 2, . . .r , belongs to Z{su{n)) , if 
and only if Q^...,^ are coefficients of the polynomial in ^1,^2, ■ ■ ■ variables 

which is invariant under the adjoint action 

0(^1, 6, • • • , ^r) = 0((Ad gf^,, (Ad gf^2, ... (Ad gf^r) , 
with (Ad g)'^- the matrix of adjoint operator. Ad g , calculated in the basis 

Cil 612 . . . ^ir ■ 

Therefore, from the algebraic standpoint, the study of the orbit space 
|SU(n) as well as any characteristics of quantum-mechanical observables, 

invariant under the unitary action (12. ip . reduces to the computation of the 

center Z(su(n)) of it(su(n)). 
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If the elements (Er belong to Z they are termed as Casimir operators. The 
number of independent homogeneous Casimir generators for SU(n) group is 
equal to rank su(n) = n — 1 . 

It is well known, that the quadratic Casimir operator is unique up to the 
constant factor and is expressible with the aid of the Cartan tensor: 

a, =tr((Ad e,)(Ad e,)) , 

Therefore for algebra 5u(n) the quadratic Casimir operator reads 



The higher dimensional Casimirs are expressed via the symmetric struc- 
ture constants dijk of su(n) algebra p!5]. Because further, dealing with the 
qubit-qutrit system the Casimirs of SU(6) will be used[§, the expressions for 
(Tj are given below: 



£6 



^ ^ dii-i^i2diji^djj^i^dki^i^ 6ij^6i2€i^6i^6i^6i^ 



Now using these operators and decomposition fl2.3p based on the iso- 
morphism between center Z{su{n)) and SU(n)-invariant polynomials, the 
following scalars, referred hereafter as Casimir invariants, can be written: 

(2.6) 
(2.7) 
(2.8) 

(2.9) 

(2.10) 

where 

(UWV)a:= KdabcUbVc, 



^2 = 


(n 




^3 = 


1 — 


i)(^vo-^ 


u = 


(n — 


i)(^vO-(|vO 


12:5 = 


(n — 


1) ((Ivov(^vo)-C 


u = 


(n — 


1) {iy^yif 



^The tensorial su(2) ® su(3) product type basis for su(6) is given in Appendix [X] 
The formulas for the symmetric structure constants dijk and the antisymmetric structure 
constants fijk for su(n) are presented there as weU. 
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with normalization constant n := A/n(n — l)/2 . 

Now, using these scalars, the positive semi-definiteness of density matrices 
for an arbitrary n-level quantum system will be formulated. 

• Positivity of density operators • To the best of our knowledge 
the first analysis of consequences of the constraints on the density operator 
due to its semi-positive definiteness has been done in the sixtieth of the 
last century studying the production and decay of resonant states in strong 
interaction processes |T6l [TTl [18] . Nowadays, the quantum computing and 
theory of quantum information reveals the new role of these constraints and 
recently they have been once more derived [HI [20] |. 

To formulate the semi-definiteness let us choose the Bloch representation 
for a density operator 02.21) [TT] : 

^=-(In + w), U = K$,-\, (2.11) 

n 

where (n^ — 1) -dimensional Bloch vector ^ G M"^^^^ is contracted with el- 
ements Aj , i = 1^. . ,n^ — 1 of the Hermitian basis of su(n) Lie algebra. 
According to [17] |j a necessary and sufficient condition for the Hermitian 
matrix to be positive is that the coefficients 5*^ of its characteristic equation 

\lx-g\=x''- Six''~^ + S2x''-^ - ... + {-lyS^ = (2.12) 

should be non-negative 

^>0 ^ Sk>0 k = l,...,n. (2.13) 

It is convenient to rewrite these inequalities in terms of normalized coeffi- 
cients Sk := S'fc/max{ 5"^} . Since the maximal values of Sk correspond to a 
maximally degenerate roots; Xi = X2 = ■ ■ ■ = x^^ = l/n of the characteristic 
equation f l2.12p one can express them via the binomial coefficients 

-^{^4 = ^(/^) 

and thus 

0<Sk<l fc = 2,...,n. (2.14) 

■^In our recent publication [5^ the positivity conditions for density operators has been 
analyzed in context of the consequences for integrity basis of SU(2)(8)SU(2) polynomial 
invariants ring as well as for entanglement characteristics of mixed qubit states [2T| . 

''Note that P.Minnaert attributed the same result to D.N.Williams. 
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Now we are ready to rewrite the constraints fl2.14p in terms of the Casimir 
invariants fl2.6l) - fl2.10l) . This is possible since, each of three sets, Cfc, or 5*^ , 
or tk = tr {g^) , k = 2, . . . ,n forms the basis of algebraically independent 
invariants of SU(n) group (see e.g. [15]). The expressions for the coefficients 
Sk in terms of tm sue well-known, they are given by determinants: 
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Further, tm can be represented as polynomials in Casimir invariants. Based 
on the expressions for traces of symmetrized products of Lie algebra basis 
elements (see Appendix [Xj cf. also [20]) we have: 



tr (uj') =nC2 , 

tr (io^) =nC3 , 

tr (u^) =n {Cl + , 

tr {co') =n (2 £2^3 + CTs) , 

tr (lo^) =n {€l + 2 ^2^4 + ^3 + U) ■ 

Finally, imposing the following normalization for the Casimir invariants. 



n 



l)(n-2) 



k + l) 



we arrive at a system of inequalities in su(6) Casimir invariants, that defines 
the positive semi-definiteness of the density matrix of qubit-qutrit pair: 

0<C2<1, (2.15) 

0<3C2-C3<1, (2.16) 

< 6C2 -5^2^ -4C3 + C4 < 1, (2.17) 

< (1 - 5 Cif - 30 C2C3 + 10 C3 - 5 C4 + C5 < 1 (2.18) 

<(1 - 5^2)^ - ISQCiC^ + 125C2C4+ 

+ 20^3(1 + 5C3) - I5C4 + 6C5 - Ce < 1 . (2.19) 



Since the positive semi-definiteness of density matrices plays an excep- 
tional role in the entanglement quantification problem it is reasonable to 
rewrite the system of inequalities (I2.15p - (l2.19p in terms of the local SU(2) ® 
SU(3) invariants. 

3 The local unitary invariants 

• The local invariance of composite states • When a quantum system 
is obtained by combining of r-subsystems with ni, n2, . . . , n^ levels each, the 
non-local properties of the composite system are in correspondence with a 
certain decomposition of the unitary operations f l2.ip . 

In order to discuss this decomposition consider the subgroup of unitary 
group formed by the local unitary transformations 

SU(ni) SU(n2) ® ■ ■ • ® S\J{n,) , (3. 1) 

acting independently on the density matrix of each subsystem 

Q^^^^ ^ 0^""^^' = g0^''^^g-' geSVin,), t = l,2,...,r. 

Two states of composite system connected by the local unitary transfor- 
mations fl3.ip have the same non-local properties. The latter can be changed 
only by the rest of the unitary actions, 

SU(n) 

SU(ni) ® SU(n2) ® ■ ■ ■ ® SUK) ' ~ ""'"^^ 

generating the class of non-local transformations. 

Having this notions, we are in position to discuss the structure of the cor- 
responding ring of polynomial local invariants, i.e. polynomials in elements 
of the density matrices, which are scalars under the adjoint local unitary 
transformations. It is well known that for any reductive linear algebraic 
group G and for any finite dimensional G-module V, the ring has the 
Cohen-Macaulay property [22] and possesses the Hironaka decomposition 

r 

a=0 



8 



where i^^^, 6 = 1, 2, . . . , s are primary, algebraically independent polynomials 
and J a, a = 0,1,2, ... ,r, Jq = 1, are secondary, linearly independent in- 
variants respectively. According to that the corresponding Molien function 
Mci^q) for TZ'^ [7] can be expressed as follows 



In this form it provides us with a certain knowledge on the numbers of alge- 
braically independent polynomials as well as linearly independent invariants. 

• Molien functions for Cl^^^^"^^] and C[*pJ®^^] • Let us start with 
a remark concerning the adjoint action fl2.ll) . Consider a non-degenerate 
density matrices. In this case using the natural identification of the elements 
of a linear space spanned by the Hermitian n x n matrices with the space 

g pij 

one can instead of the adjoint action ( 12. ip consider the linear representation 
on M°'-i 



= LabVb , Lab e SU(n) ® SU(n) , 

where the line over expression means the complex conjugation. 

After this identification in order to get some insight on the structure of 
the ring of polynomial invariants of linear action of Lie group G on the linear 
V space we can compute the Molien function 

M{C[Vf, q)= [ ,y |g| < 1 , (3.2) 

Jg det(I - qn{g)) 

where dfi{g) is the Haar measure for Lie group G and TT{g) is the correspond- 
ing representation on V . We start with the system of two qubits. 
Two qubits . In this case the local unitary group is 

G = SU(2) ® SU(2) . 

As it is well known for any reductive linear group the integration in fl3.2p 
reduces to the integration over the maximal compact subgroup K of G [1]. 
In the present case this results in integration over the maximal torus 



■K{g) = diag {1,1, z, z ^) ® diag {1,1, w,w ^) 
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where z,w - coordinates on one- dimensional tori. Therefore computations 
reduce to the following two-dimensional integral 

where 

det{I -qn{g)) = {1 ~ q) ^{z,w,q) 
^(z, w, q) =(1 - qf{l - qzf{l - qwf{l - qz-^f{l - qw'^f 
(1 — qzw){l — qz~^w){l — qzw^^){l — qz~^w~~^) . 
After integration we get the Molien function [7] 

1 + + + 3g6 + + 2q^ + Sg^ + + q^^ + q^^ 

Msu(2)^SU(2)(g) - (1 -g2)3(l_ ^3)2(1 _ ^4)3(1 _ ^6) ' 

which is the palindromic one, 

^SU(2)iX)SU(2)(l/q') = -Q'^^^SU(2)(g)SU(2)(^) , 

with the degree consistent with 

dim SU(4) = 15 . 

Qubit-Qutrit . Now the local unitary group is G = SU(2) ® SU(3) and 
owing to the symmetries of the integrand (13. 2p the non-trivial part of the 
integration is entirely accumulated in the diagonal components of the T^{g)- 
representation of the form: 

7r(5')diag = diag (1, 1, x, x~^) ® diag (1, 1, 1, y, z, yz, z~'^, iyz)~'^) , 

where x, y and z are coordinates on the maximal torus. Therefore, the com- 
putation of the Molien function (13. 2p reduces to the evaluation of the multiple 
contour integral in complex planes over unit circles 0: 



Msu(2)®su(3)(g) = TTT-T^ / / / f{x,y,z,q)dxdydz , (3.3) 
(27rz)'^ 



'^The multiple integral p.3p has been calculated using the consecutive application of 
Cauchy's residue theorem. Since the integrand / has poles of rather high orders, com- 
puter computations of the residues has been performed using the command Residue built- 
in Mathematica that implements the standard limit formula for high order poles (see 



http : //mathworld . wolfram . com/ComplexResidue . html ). 
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where 

xyz '^{x,y,z,q) 
det(I - qn{g)) = {1 - q) ^(x, y, z, q) , 

and y, z, g) = 

(1 - qf{l - - g.)^(l - qyzf{l - ^)^(1 - ^)^(1 - 

(1 — qx^il — qxy){l — qxz){l — qxyz){l — — )(1 — — )(1 — — ' 

1/ z yz' 

(1 _ 1)3(1 _ _ _ _ _ ^ ^ 



a; a; a; x xy xz xyz 

As a result, the Mohen function can be represented in the rational form 

(cf. my- 

N 

^SU(2)®SU(3)(g) = , (3.4) 

where 

AT =1 + 4 + 9 + 38 + 69 + 173 g^ + 347 g^ + 733 g^" + 1403 g^^ 
+ 2796 g^2 + 5091 g^^ + 9286 g^^ + 16058 g^^ + 27208 g^^ + 44250 g^^ 
+ 70537 gi^ + 108430 g^^ + 163158 g^^ + 238264 g^^ + 339974 g^^ 
+ 472130 g^^ + 641187 g^^ + 848615 g^^ + 1098643 g^^ + 1388741 g^^ 
+ 1717327 g^® + 2075836 g^^ + 2456389 g^° + 2843020 g^^ + 3222408 g^^ 
+ 3575226 g^^ + 3884797 g^^ + 4133599 g^^ + 4308636 g^^ + 4398377 g^^ 
+ 4398377 g^^ + . . . + 38 g^^ + 9 g™ + 4 g^^ + g^^ 

D =(1 - g2)3(l - g3)4(l _ g4)5(^ _ ^5)4(^ _ ^6)5(^ _ ^7)2(^ _ ^8) _ 

This Molien function is the palindromic one 

A^'sU(2)®SU(3)(1/Q') = Q'^^^SU(2)®SU(3)('?) , 

as provided by 

dim SU(6) = 35 . 

This form of the Molien function serves as source of information on the 
polynomial ring of SU(2)(8> SU(3) invariants. Particularly, one can endeavour 
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to identify the structure of algebraically independent local unitary scalars. 
According to fl3.4p there are 24 independent scalars in agreement with simple 
count of dim[SU(6)/SU(2) SU(3)] = 35 - 11 = 24. The set of these 24 
polynomial invariants may be composed from three invariants of degree 2, 
four of degree 3, five of degree 4, four of degree 5, five of degree 6, two of 
degree 7 and one of the degree 8. 

Note that the Poincare series of Msu(2)®su(3)(?) 

oo 
(1=0 

determines the number of homogeneous polynomial invariants of degree d. 
According to the calculations of fl3.3p the few terms of the Taylor expansion 
over q are 

Msu(2)®su(3)(g) = l + 3g2 + 4g3 + i5g4 + 25g^ + 90g^ + 170g^ + 489g« 
+ 1059 + 2600 + 5641 q^^ + 12872 q^^ + 27099 q^^ 
+ 57990g^^ + 118254 + 240187 + (g^^) . (3.5) 

Now, having in mind the input from the structure of the Molien function 
(13. 4p . we attempt to construct the local SU(2)®SU(3) unitary invariants. 
• Constructing SU(2)(8)SU(3) invariants • Let us introduce the de- 
composition for density matrices well adapted to the case of composite qubit- 
qutrit system. The space su(6) in (12. 2p for n=6 admits decomposition in the 
direct sum of three real spaces 

3 

su(6) = K = 5u(2) ®h + h® su(2) + su(2) ® su(3) . 

a=l 

Using Pauli matrices cTj as the basis for su(2) and Gell-Mann matrices as 
the basis for su(3) (see Appendix|X]) the density matrix (12. lip for qubit-qutrit 
system can be written as [9l [TO] : 

g=^ [Iq + u] , w = a + /5 + 7, 
6 

where 

3 8 3 8 

a := ^aiO-j® I3, /3 := ^bah K , 7 ■= ai \a . 

i=l a=l i=l a=l 
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Among the 35=3+8+24 real parameters (oj , ba, Cia) the first two sets, Oj 
and ba , correspond to the Bloch vectors of an individual qubit and qutrit 
respectively; the evaluation of partial trace yields the reduced matrices for 
subsystems: 

■.= tiB{g) = ^{h + a-a), ■.= tTA{g) = ^ih + b-X), 

while the variables Cia are entries of the so-called correlation matrix C = (cia). 

Now using a trace operation described below we can construct a set of 
local SU(2)(g)SU(3) scalars, candidates for the elements of the integrity basis. 

In analogy with the generators (12 .Sp of the universal enveloping algebra 
we consider a set M of non-commutative monomials 

•= ■ ■ ■ ■ ' G Ai , 

where each of Xi^ ,k = 1, . . . ,d , is one of a , /3 , or 7 . To each we 
assign a multidegree {s,t,q), s + t + q = d , where s, t and q are degrees of 
a , /3 , and 7 respectively. The trace operation on monomials 

tr : Mi^...i^ Vstqicii ,ba, Cia) ■■= ti {Mi^...iJ e 'P ■ 

defines the map tr : M ^ V of into the algebra V of homoge- 
neous polynomials in variables (aj , ba, Cia)- A generic term of the polynomial 
'Pstqi^'i , ba, Cia) IS a couvolutiou of vcctors ttj , ba and matrix Cia with traces 

tr {aia2 ■ ■ - (Tp® A1A2 ■ ■ ■ A,.) = tr (crio-2 ■ ■ - c^p) (A1A2 ■ ■ ■ A^) , 

where p = s + q and r = t + q. 

Now it is easy to verify that the images of the trace map are local unitary 
invariants. Indeed, since under the transformation of the form fci ® ^2 , where 
ki G SU(2) , and ^2 ^ SU(3) , the matrices ex's and A's in the basis elements 
of 5u(6) (see Appendix |A]) are transformed independently, in adjoint manner 

a — )■ kiaki^ , A — t- k2\k2^ , 

the polynomials tr {Ai) are invariant against SU(2)(8>SU(3) action. 

Therefore the polynomials Vstqi^i , ba, Cia) are the reserve for the integrity 
basis of the ring C[^]^^(^)'^^^(^) . Now, in contrast to the case of SU(n) 
Casimir invariants built up with the help of symmetric structure constants 
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only, dealing with the scalars against the tensor product of groups, the in- 
variants are constructed in terms of the antisymmetric structure constants 
as well. For example. 



tr(7^) = QaCjfeCfectr {aiGjak ® XahK) = CiaCjbCkctr {(Jiajak) tr [XaXbK) ■ 



This quantity being invariant under the SU(2)(8>SU(3) action is expressible 
via totaly antisymmetric tensor eijk - structure constants of su(2) algebra 
and fabc - structure constants of su(3): 



Choosing a basis for local invariants, several types of algebraic depen- 
dence between the polynomials in V have to be taken into account. It is 
worth to consider two illustrative examples. Applying the Hamilton-Cayley 
theorem for elements a, P and 7, considered as Hermitian 6x6 matrices, 
one can determine the algebraic identities for the polynomials of the degree 
d > 7 . Less obvious example of relations between polynomials is due to 
the identities between the structure constants of the algebra. |^ Let us con- 
sider two invariants, both 4-th order in variables of the correlation matrix 
C , but one constructed using the invariant symmetric structure constants d 
while the second one with the anti-symmetric structure constants / : 



With the aid of identities flA.ip and flA.2p (Appendix for the structure 
constants of su(3) algebra, one can convinced that 



According to the Poincare series (13. 5p there are 15 homogeneous scalars 
of order 4, while there are 81 = 3^ monomials in three noncommutative 
variables. But since the elements a and (3 commute this number reduces. 
Taking into account this commutativity as well as the invariance of trace 
operation under the cyclic permutations of products, one can find 18 valuable 

^For the detailed analysis of the relations of that type we refer to [23] . 



tr(7 ) — ^^ijkfabcCiaCjbCkc- 



3^'^\dd) 




3°^\dd) = - J^'^^ff) - - [(tr(C^C))' - 2tr(C^CC^C) 
3 3 - 
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monomials: 

a^ /3^ 7^ a^^ , afi^ , a^'j , 07^ , /3^j , 13'^^ , 

a^/37 , a/9^7 , a/37^ , q;7/37 . 

Taking traces of these monomials one can convince that five of them form 
the kernel of trace map: 

tr(a^/3) = ti{a(3^) = tr{a^-f) = tr{(3^-f) = tr{a^l3-f) = , 

and images of last two monomials coincide up to sign 

tr (a/37^) = — tr (07/37) . 

Therefore the following set of twelve traces 

tr(a^), tr(/34), tr (a^/?^) , tr (a V) , 
tr(7^), tr(a7=^), tr (/37=^) , tr (a7a7) , 
tr (Z^^^) ^ tr (/37/37) , tr , tr (a/37') , 

plus three 4-th order polynomials constructed as products of second degrees 
polynomials tr (a^)tr (/3^) , tr (a^)tr (7^) , tr (/3^)tr (7^) , are 15 homogeneous 
invariant polynomials in accordance with the Poincare series (13. 5p . 

How difficult is it to extract the independent scalars from this list? It is 
easy to verify that some traces are expressed in terms polynomials of second 
order; e.g., tr (a^/3^) = |tr (a^)tr (/3^) . Concerning the remaining monomials 
one can see that several of them have the same multidegree. Namely, the 
following "trace" polynomials 

1. tr (q;^7^) = |tr (a^)tr (7^) and tr (Q;7a7) , 

2. tr (/32)tr (7^) , tr {(3^^^) and tr {(3^P^) , 

belong to the same space V202 and 7^022 respectively. Being linearly indepen- 
dent monomials, they obey the following relations 

tr (a^7^) + tr (0707) = 8 ai.ai^Ci^j^Ci^j^ , 

tr {f3^-f^) - itr (/3^)tr (7^) = 4 dj^j^k dkj.j, bj.bj^d.j.d.j^ , 

tr (/3 7 ) + tr ((3-y(3'y) = 8{-^bj-^bj^Ci-^j-^Ci-^j^ + dj-^j^^. d^j^j^ bjJjj^Ci-^^j^Ci-^j^ , 
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where summation over all indices is assumed. This circumstance leaves an 
open question how to build the elements of integrity basis with a certain 
multidegree using the "trace" polynomials. 

We resume our analysis by the following list of linearly independent 
SU(2)®SU(3) scalars which are not products of low orders ones0: 

• degree 2, three invariants 

tT{a'), tr(/32), tr(72), 

• degree 3, four invariants 

tr(/33), tr(73), tr(a/37), tT{/3j') , 



• degree 4, eight invariants 

tr (7^) , tr (073) , tr {(3^^) , tr (a7a7) , 

tr {(3^^^) , tr (/37/37) , tr (a/?^) , tr {a^^') . 

• Casimir invariants decomposition • The expansion of the SU(6) 
Casimir invariants up to the 4-th order fl2.6p - fl2.8p over the above suggested 
SU(2)®SU(3) "trace" scalars reads: 

6 €2 = tr (a^) + tr + tr (7^) , 

6 €3 = tr {/S"^) + tr (7=^) + 3 tr {(3^^) + 6 tr (a/37) , 



1 

3 

+ 4 



tT{a') (2 tr(/32) + tr(7^) j + \ tri^^ - i tri^^ - triP'M^') 
tr(a7^) + tr(/373) + tr(/3 V) + tT{a(3-f^) + 3tr(a/3%) 
tr(a7a7) + ti((3'j(3'j) + tr(7^) . 



+ 2 

We conclude with the final remark on the applicability of the derived re- 
sults to the problem of classification of mixed quantum states. Using inequal- 
ities f|2.15p - fl2.19p and results from [21] the well-known Peres — Horodecki cri- 
terion for the separability of qubit-qutrit mixed states can be reformulated 
as a set of inequalities in SU(2)®SU(3) scalars. 



^Note that 2-nd and 3-d order invariants were proposed in 



16 



• Acknowledgments • The work was supported in part by the RFBR 
(grant No. 10-01-00200), and by Ministry of Education and Science of the 
Russian Federation (grant No. 3810.2010.2) and JINR-Bulgaria 2010 collab- 
orative grant. The research of A. K. was supported by the GNSF research 
grant GNSF/ST08/4-405. 

A Appendix: Formulas for the 5u(6) algebra 

• The tensorial basis • For the su(6) algebra we use the basis {ta}a=i,...,35 
constructed from the tensor products of the Pauli matrices Uj e 5u(2): 

01\ f -i\ /ID 



and eight {Aa}a=i^...^8 Gell-Mann matrices, forming the 5u(3) basis: 
Ai = 



A4 




The elements ta are enumerated as 

1 . 1 . , 

= (Tj (8) I3 , Ts+a = -^h<^K, 

Tll+a = ^ (T"! ® Aa , Tig+a = (72 <H) Aa , T27+0 = ^ CTS <8) Ao . 

• The algebraic structures • The product of basis elements reads 

2 

taTb = -Sab'^ + {dABC + I fABc)rc , 
n 
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The structure constants dABC and Jabc can be determined via equations 

1 2 
dABC = -Tt{{ta, tb}tc), Iabc = --Tr([rA, rB]rc), 

where apart from the Lie algebra product, [ , ] , the "anti-commutator" of 
elements, i.e., {ta, tb} = taTb + tbTa has been used. 

• Identities for structure constants • For the SU(n) group the the 
structure constants obey the following identities: 

fabcfcpq ~l~ fbpcfcaq ~l~ fpacfcbq ; 

2 

fabcfcpq dapcdcbq duq^dQ^jp -\- — {^ap^bq ^aq^bp) i 
fabcfcpq ~l~ faqcfcpb '^dg^pQdf.ffq d^bcdcpq '^agc'-^cftp 

2 , 

+ -{26apSbq - SabSpg - Saq6bp) . (A.l) 

The SU(3) symmetric constants satisfy [211 EH] an important identities 



dabcdcpq ~\~ d^pcdcaq ~\~ dpacdcbq g i^^ab^pq ~\~ ^ap^bq ~\~ ^aq^bp) ■ (-A'^) 

• The traces • The traces of symmetrized products of su(n) basis 
elements are 

tr ir^aU}) =2 Sab, 

tr {r^aTbTc}) = 2dabc, 

2^ 

tT {TJaTbTcTd}) = Sabred + 2 dabedecd , 

n 

2^ 

tr {r^aTbTcTdTe}) = — {dabc^de + ^abdcde) + 2 dabfdfcgdgde , 
23 2^ 

tr {T{aTbTcTdTeTf}) = — Sab^cd^ef + — {dabgdgcdSef + ^abdcdgdgef) + 



~l~ ^ dabcddef ~l~ 2 dabgdgchdhdvdvef 
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